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Stress-Strain Relationships in Yarns Subjected to Rapid 

Impact Loading: 5. Wave Propagation in Long 

Textile Yarns Impacted Transversely 

Jack C. Smith, Frank L. McCrackin, and Herbert F. Schiefer 

The behavior of an infinitely long flexible filament after transverse impact is treated 
theoretically. The filament is assumed to have a tension-strain curve that is always con- 
cave downward, and to have no short-time creep or stress-relaxation effects. Under most 
conditions the impact initiates a variable strain that propagates down the filament between 
an "elastic wave" front and a "plastic wave" front. A transverse wave, shaped like an 
inverted V, then travels in the constant-strain region behind the plastic-wave front. Under 
special conditions the transverse- wave front may propagate faster than the plastic- wave 
front, but the shape of the transverse wave remains the same. The theory for both cases is 
worked out in detail, and some illustrative examples are given. 

1. Introduction 

The behavior of textile varus held in chimps 40 to 60 cm apart and impacted transversely 
was demonstrated and discussed in the fourth paper of this series [1, 2, 3, 4]. 1 The present 
paper deals with the theory of transverse-wave propagation in yarns of infinite length where 
there are no clamps to reflect the waves. The theory developed here will be used in future dis- 
cussions of transverse impact. 

When a yarn is struck transversely, a strain wave is propagated along tin 4 yarn outward 
in each direction from the point of impact. In the region between each of these wave fronts, 
material of the 4 yarn is set into motion longitudinally toward the point of impact. This ma- 
terial is taken up by a tent-shaped wave of transverse motion that propagates outward from the 
impact point. This behavior is demonstrated in figure 1, in which several configurations of a 
yarn after transverse impact are shown. The position of the longitudinal-strain wave in each 
of the configurations is shown by the arrows. The backward flow of yarn material is indicated 
by the white tick marks. 

/^'^^\\\ Figure 1. Configurations of a yarn after transverse impact. 

^/O^^^On^^v, These configurations, 40X10" 6 sec apart, were calculated for a yarn having a 

/ / /S N \X X N \ \ X N linear (Hooke's law) tension-strain curve. Longitudinal-strain wave velocity 

^/^ N \ \V\ V \ N ~~~~~~~~ ~~ was 2 ' tW0 m/sec. The positions of the longitudinal-wave front in successive con- 

"'VO^^O-o^ x \- figurations, are designated by the arrows. Impact occurred at the vertex of the 

/ /'^\^\\^\^ N _ri~~~~~~~~~~3~~~~~~~~~~~" transverse wave shown at the left, at impact velocity of ISO m/sec. The leading 

^^^^\\. \^\^— — ~~~~~~~~~Z1Z:~~~~~~~~~'Z:~~ C(1 S e of tne transverse wave propagates with velocity '-M) m/sec with respect to 

""V^'O^-s^/*' ZL~~~~~~~~~~~~~Z:Z1~~'Z:~~~~ the uns trained yarn material. Inward flow of yarn material at velocity 50 m/sec 

/ /^>\\ Vv \ J is shown by the white tick marks. There is no horizontal component of flow in tho 

/ /'^\^^^~~~~~~~~~~~~~~~~~~~~ J Z:~~~Z*'ZZZ:Z region of the transverse wave. 
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The configurations shown in figure 1 were calculated for a yarn composed of a Hookean 
material, using relations developed in this paper. The longitudinal strain-wave velocity was 
2,500 m/sec. Transverse-impact velocity was 180 m/sec. Under these conditions a 2-percent 
strain is produced in the wake of the longitudinal-wave front, and yarn material in the region 
between longitudinal- and transverse-wave fronts flows inward at a velocity of 50 m/sec. In 
the wake of the 1 ra nsverse-wave front, which propagates at a velocity of 350 m/sec, the yarn 
material flows transversely at the impact velocity as shown. The time interval between each 
configuration is 40X10" 6 sec. 

Transverse waves of the type just described have been discussed by Taylor and his associ- 
ates [5]. Cole, Dougherty, and Huth [6] have treated the problem of transverse-wave propaga- 
tion in prestressed wires. Craggs [7] and McCrackin [8] have also discussed the problem for 

1 Figures in brackets indicate t tic literature references at the end of this paper. 
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plastic-elastic strings, showing the effect of wave reflections and interactions. The problem 
of oblique impact of an elastic cable has been discussed recently by Ringleb [9]. The treat- 
ment presented here extends the work of these authors and makes use of concepts introduced 
by von Karman [10] and Taylor [11] in their treatments of longitudinal impact. 

2. Theory of Longitudinal Impact 

Longitudinal-strain waves resulting from a constant-velocity tensile impact do not neces- 
sarily have sharp wave fronts. According to von Karman, the strain increases gradually in 
the region between two wave fronts, an "elastic wave" front 2 propagating at velocity 



^=^=Vm(S1o (i) 

and a "plastic wave" front 2 propagating more slowly at velocity 



^=^)=V/X£L; <*> 



In eq (1) and (2), M is the mass per unit length of the unstrained filament material, and dT/de 
is the slope of the tension-strain curve, evaluated at zero strain in eq (1) and at strain e p in 
eq (2). In the wake of the plastic-wave front the strain in the filament has the constant 
value t p . The value of e p depends upon the speed at which the filament is impacted longi- 
tudinally. 

The filament material in the wake of the plastic-wave front flows in the impact direction 
at a velocity uniform along the filament given by 



H 4VB(f>=J> w '- (3 > 

For a filament impacted longitudinally at one end, W is the same as the impact velocity. 

In deriving the above equations, von Karman assumed a tension-strain curve that was 
always concave downward. Tension-strain curves for textile yarns often have a concave up- 
ward portion which, however, is less pronounced in curves obtained at rapid testing rates. In 
addition, textile yarns exhibit creep and stress-relaxation phenomena. Such phenomena were 
not considered by von Karman. Despite this, von Karman's treatment represents a notable 
advance over theories based on Hookean tension-strain behavior. 

3. Theory of Transverse Impact 

In the laboratory or observer's coordinate system, let x' and y' denote the horizontal and 
vertical coordinates, respectively. Let a filament lying along the x' axis be impacted with 
velocity V in the y' direction at the point where x' = 0. The motion of the filament is to be 
described in terms of Lagrangian coordinates. That is, distances along the unstrained length 
of the filament are given in an x coordinate system, which before impact coincides with the 
x f coordinates. The displacement of any point x of the filament at any time t relative to its 
original position is given by £(x,t) for horizontal displacements and rj(x,t) for vertical displace- 
ments. The location in the laboratory coordinate system of a point x on the filament is thus 
given by 

x'=x-\-£(x,t), (4a) 

y'=n(x,t). (4b) 



2 The terms ''elastic" and "plastic" are used as convenient designations for the wave fronts bounding the region of variable strain. They are 
not intended to characterize the complex recovery behavior of the material. For example, the theory holds for materials having nonlinear tension- 
strain curves, but which recover almost completely from large deformations. For these materials, the strain associated with the plastic wave is 
almost entirely elastic. 
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Consider an element of the filament of unstrained length Ax. After impact, this element is 
strained, and its new length, measured in the laboratory coordinate system, becomes 
Aa:V [1+ (d£/dx)] 2 + (d^/dx) 2 . The increase in length per unit length, or strain e, is thus 



•=V( i+ 0+@)*- 1 - 

The~angle that the element makes with the horizontal is given by 

™ e =wM> < 6a) 

cos. = ^(1+1) (6b) 

The differential equations that govern the motion of the smalTelement are 

*$-&?«*<>' (7b) 

where M is the mass per unit length of the unstrained filament, andJTis the tension. By making 
use of (6a) and (6b), these equations become 

vt,_br t an 

Assume, as von Karman did, that the tension Tis a function T(e) of the strain only, given 
by the tension-strain curve. For solutions in which the strain, and consequently the tension, 
are constant along the filament, eq (8a) and (8b) reduce to 

d^_ T d 2 7? . 

W~M(l + e)c)x*' {y&) 

W M(l + e)dx*' {V ' 

£ and 7) must also satisfy eq (5). 

For other solutions in which ij=0 (no transverse wave), eq (8a), (8b), and (5) reduce to 

t?=0, (10a) 

W = M*-x' < 10b) 

This is the case considered by von Karman, and to which his solution, eq (1), (2), and (3), is 
applicable. 

A combination of solutions for the preceding special cases gives the solution to the problem 
in which a long filament initially at rest along an axis in the x' direction is impacted transversely 
at x=0 with velocity V. The boundary conditions that must be satisfied for the x^>0 half of 
the filament are £=0 and rj=Vt at a;=0, and £=r}=0 at x= co. 
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4. A Simple Solution Usually Applicable 

Immediately after impact a train of longitudinal-strain waves propagates outward along 
the filament. The leading wavelet in the train, or elastic wave, propagates at velocity C e given 
by eq (1). At points in advance of the elastic wave the filament is unstrained. The elastic 
wave is succeeded by a series of other more slowly propagating wavelets, each wavelet adding 
an increment to the strain in the filament. In the wake of each wavelet, material of the fila- 
ment flows inward toward the point of impact at a velocity that increases as each wavelet 
passes. The. velocity of each wavelet is proportional to the square root of the slope of the 
tension-strain curve, hence in order that the wavelets do not overtake each other, the slope 
of the tension-strain curve must decrease as e increases (be concave downward). The final 
wavelet in the train, or plastic wave, propagates with velocity C p given by eq (2). The strain 
in the wake of the plastic-wave front has the constant value e p . Filament material in the wake 
of the plastic-wave front flows inward at the constant velocity W given by eq (3). 

For the simple solution discussed here, transverse-wave motion of the filament occurs only 
in the wake of the plastic-wave front. The transverse wave thus propagates with velocity Z7, 
which is less than C p . For the region in advance of the point Ut on the filament, where only 
longitudinal-strain waves occur, von Karman's solution [10] is applicable. This solution, 
satisfying the differential eq (10) and the boundary conditions £= — Wt at x=Ut and £=0 
at x= co } may be written as follows: 

For the region O e t<Cx<C °° , 

-f=*=6=0. (11) 

For the region C p t<x<C e t, 

77=0. (12a) 

To find the amount of motion in the £ direction, consider a wavelet, which increases the 
strain in the filament by an increment de. This wavelet propagates at velocity C(e) and arrives 
at point x at time t=x/C(e). The inward velocity of material flow is increased by an amount 
— C(e)de during time t-xlC(e). The distance £ that point x moves inward is the integrated 
product of the velocity and time after arrival for each wavelet passing through the point x, or 



where 



t=-£0(a)[t-^^da=-W(e)t+ex, (12b) 



*<WaKS> (l2c) 

W(e)==( § C(a(da), (12d) 



and a is a variable designating strains along the filament in advance of the point where the 
strain is e. 

During the time interval between the arrivals of the elastic and plastic waves at point x, 
the strain is found as a function of x and t by solving the following equation for e: 

G(e)=x/t. (12e) 

For the region Ut<x< C p t, 

71=0, r (13a) 

e=€p=a constant. (13c) 

In order to describe the motion m the transverse-wave region (0<x< Ut) it is necessary to 
find solutions of eq (9a), (9b), and " (5),. which satisfy the appropriate boundary conditions.. 
The most general solution of (9a) ha s the form ri=f+g-{- At+Bx+K } where/ and g are non- 
linear functions of the variables t—x/^T/M(l + e) and t+x/-yjT/M(l + €); respectively. Applying 
the boundary condition that i^Wai x=0, gives Vt=j(0,i) + g(0,t)+At-\-K. In order that 
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this relation be satisfied, it is necessary that K=0 and g(0,t) = - 
must be the same function except for sign, and 



-f(0,t). Therefore,/ and g 



v=At-xMT/M(l + e))-f(t+xMTIM(\ + e)) + Vt+Bx. 

The condition that rj=0 at x=Ut can only be satisfied if /=0 and B= — V/U. Thus 
rj=V(t—x/U). For this value of 77 and for a constant value of e, the most general solution of eq (5) 
becomes £=:z[V(l + e) 2 — (V/U 2 — l]+h(t)+K. The condition {=0 at x=0 requires that 
h(t)+K=0. The solution for £ just obtained is also a solution of eq (9b); thus: For the region 
<x<Ut, 

v =V(t-x/U), (14a) 



i:= x y(l + e 2 ¥-(V/U) 2 -l], 
€=e 2 =a constant. 



(14b) 
(14c) 



The behavior of the negative half of the filament is the same as that just described except 
for the direction of wave propagation. 

The configuration of the positive half of the filament at time t after impact is shown in 
figure 2 by the line passing through the points A, B, P, Q. Filament material in the transverse- 
wave region 0<x< Ut, moves in the y r direction because d£/d£=0 in this region. Thus points 
A, B, P moving with velocity V, arrive at positions A*, S*, P* at the time t-\-dt. The transverse- 
wave front moves with velocity U relative to points on the unstrained filament. The horizontal 
component of the velocity relative to the x f , y f coordinate system therefore is Z7(l + e 2 )cos 0, 
and the wave front at point P advances a distance Udt(l-\-e 2 )cos 6 to point Q* in a time dt. 
Filament material in advance of the transverse-wave front is moving inward with velocity W, 
so that in time dt point Q moves a distance — Wdt to its new position Q*. 




— J wdt U- 



Figure 2. Configuration of positive half of filament l seconds 
after impact. 

Filament originally lay along the x' axis and was struck with velocity V in the 
y' direction. After t seconds the positive half of the filament has the configuration 
ABPQ, after t-\-dt seconds the configuration A*B*P* Q*. The filament is assumed 
to have strain 62 and tension 1\ in the transverse-wave region ABP— A*B*Q*, 
and strain e p and tension T P in the region in advance of the transverse wave. (It 
is proved in the text that in most cases e2= c P , Tz= T p .) Filament material in the 
transverse wave moves vertically with velocity V; material in advance of the 
transverse wave moves inward horizontally with velocity W. The transverse- 
wave front moves with velocity U relative to points on the unstrained filament 
but the horizontal component of its velocity relative to there', y' coordinate system 
is C7(l+c2) cos 6. The point Q on the filament moves inward to position Q*, the 
location of the transverse-wave front at time t-\-dt. 



udt(i + e 2 )cos# 



The angle 6 that the transverse wave makes with the horizontal is constant along the 
filament because d£/dx and i>rj/dx are constant in the region 0<x<Ut. The transverse \\ ave is 
therefore a straight line. 6 is given by either of the relations 



sin = - 



1 V 
l + e 2 U' 



cos 



* = itW 1+£2)2 -® 2 



(15a) 
(15b) 



found by substituting tbe solutions (14) into the relations (6). 
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At the point x=Ut where solutions (13) and (14) join, £ and t\ must be continuous. Con- 
tinuity of rj is satisfied because rj=Q for both solutions. Continuity of £ requires that 



UW(l+e 2 )*-(VIUy- 1]= tJCJ-W. 



(16) 



Horizontal and vertical components of the force must also be continuous. Equating the 
horizontal components gives 



l\ 



T P -MUW=T 2 cos 6 = -^- V(l+e 2 ) 2 -(TO 2 - 



(17) 



MU is the element of mass over which the transverse-wave front passes in unit time, so MUW 
is a force equal to the change of longitudinal momentum per unit time. 
Equating the vertical components gives 



MUV=-T 2 smB = 



T 2 V 



l + e 2 U 
The velocity of propagation of the transverse wave U is found from eq (18) to be 



u= Ji 



(18) 



(19) 



'M(l+e 2 ) 
By appropriate manipulation of eq (16), (17), and (18) it is possible to obtain the relation 

The meaning of this relationship is illustrated in figure 3, where a typical tension-strain curve is 




Figure 3. Tension-strain curve whose slope decreases to a value 
less than T/(l + e). 

plotted. The long line extending from the point — 1 on the strain axis up to and through any 
point e,T on the tension-strain curve has the slope 27(1 + e). The only values of T and e 
satisfying eq (20) must be given by intersections of this line with the tension-strain curve. 
In most cases of practical interest only one intersection exists, and the relation (20) can then 
be satisfied by 



T p (e p )=T 2 (€2). 



(21a) 
(21b) 



Suppose, however, that the tension-strain curve bends over as shown in figure 3 until its 
slope for large values of e diminishes to or becomes less than T/(l + e). A second intersection 
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is possible and two sets of e,T values would be obtained satisfying (20). If a material having 
such a tension-strain curve were impacted transversely at low speeds, small values of e and T in 
the wake of the plastic-wave front would be obtained. The slope of the tension-strain curve 
for the se low values of e is greater than 27(1 + e), so the plastic-wave front would travel faster 
than -yjT p /M(l + e P ), the velocity U of the transverse-wave front. At a higher impact velocity 
a strain ei is obtained at which the slope of the tension-strain curve is just equal to Ti/(l-\-ei). 
In this special case the plastic-wave front and the transverse-wave front are at the same position 
and move together at the same velocity. If the material is impacted at a higher velocity than 
this, the solution just given no longer applies. Thus a criterion for the applicability of the 
simple solution is 

That is, the solution applies if the slope of the tension-strain curve for the material is not less 
than 7/(1 + e) for all strains less than the maximum strain e p , obtained as a result of the impact. 
Within this restriction the relations (21) can be used to obtain the results given in the following 
section. A more general solution, which includes the simple solution as a special case, is given 
later in the text. 

5. Recapitulation of the Simple Solution 

If the tension-strain curve for the filament material obtained at rapid strain rates is always 
concave downward, if the material does not exhibit large stress relaxation or creep effects 
during a small time interval, and if the conditions of the impact are such that relation (22) 
applies, the behavior of a flexible filament after transverse impact is described as follows: 

Immediately after transverse impact, longitudinal-strain wavelets propagate outward 
along the filament. The outermost wavelet, called the elastic wave by von Karman, propagates 
at velocity 



The innermost, called the plastic wave by von Karman, propagates at velocity 



Wi?(f)../ < 2 > 

P 

In the region between the elastic- and the plastic-wave fronts the strain increases to €„, the 
maximum value obtained as a result of the impact. The strain remains constant at this value 
in the wake of the plastic-wave front. 

As the strain wavelets pass a given point on the filament, material of the filament flows 
inward toward the impact point at a velocity that increases to the value 



^-/.Vra*- - (3) 

the flow velocity for material in the wake of the plastic-wave front. 

The material in the wake of the plastic-wave front forms itself into a transverse wave, 
shaped like a tent with the impact point at the vertex. The base of the tent spreads outward 
with velocity 



It should be noted that the velocity U is with respect to a Lagrangian coordinate system; 
i. e. ; a system fixed to the filament, moving and extending with it. Inasmuch as filament mate- 
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rail in the wake of the plastic-wave front is in a strained state and moving, the velocity of the 
transverse-wave front along the filament as seen by an observer in the laboratory would be 
(l-{-€ P )U— W. Similarly, the plastic-wave-front velocity in laboratory coordinates would be 
(l-\-€p)C p — W. The elastic-wave front propagates in an unstrained motionless medium and 
so would have the same velocity in both coordinate systems. 

Material in the transverse-wave tent does not move in the horizontal direction but only 
in the direction of the transverse impact. The velocity V of this motion can be expressed in 
terms of e p , U, and W by rearranging (16) to obtain 



v=Tj{i+6 P yu>-[(i+e p )u-w]> 



(24) 



If the tension-strain curve for the material is known, it is possible to calculate tables of 
values for U, V, and W as functions of e p . Then, for instance, if the impact velocity is known, 
it is possible to find the resulting strain e p and thus find the corresponding values of U and W. 

6. General Solution 

When the impact velocity is very great, the transverse-wave front in some cases pro- 
pagates more rapidly than the plastic-wave front. A solution is then required which provides 
for a variable strain in the transverse-wave region. 

Consider a tension-strain curve that is always concave downward and whose slope eventu- 
ally becomes less than Ti/(1 + €i) , where ei, T x is the point on the curve such that dT/de= 27(1 -fe) . 
Such a curve is shown in figure 4. In the region where €!<€<€*» let this curve be approximated 
by a number of straight-line segments drawn between the points ei,Tij e 2 ,2 r7 2; • . . e p ,T p . 
Point € P ,T P represents the maximum strain and tension which are attained in the material 
as the result of an impact. 

For this approximate tension-strain curve, consider a solution such as that depicted in 
figure 5, in which the transverse wave is composed of a number of straight-line segments. The 
wave front of the leading segment travels at speed U x and the strain in the segment is e^ Wave 
fronts of succeeding segments travel at decreasing velocities U 2) E7 3 , . . . U p , and the succes- 
sively increasing strains in these segments are e 2 , e 3 , . . . e p . In the region in advance of the 
transverse wave von Karman's solution applies. Thus: 



For the region C e t<ix<C °° , 



9=*=6=0. 



(25) 




€, e 2 € 3 STRAIN € 6 p 

Figure 4. Tension-strain curve assumed in derivation of the 
general solution. 
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Figure 5. 



For the region Uit<x<C € t, 



Shape of transverse wave assumed in derivation of 
the general solution. 



*?o = 0, 

£0= ~W(e)t+eX. 



(26a) 
(26b) 



In the region of the transverse wave the vertical di splacement i\ n in the nt h segment w ill 
have a general solution of the form Vn=fn(t-x/^T/M(l + e n ))+g n (t+x^TIM(l + e n )) + 
Vnit—x/U^+Aj+Kn. The boundary conditions to be satisfied are r) p =Vt at x=0, r} n =r] n -i 
at x=U n t and 771=0 at x=Uit. These boundary conditions require that the K n =0 and that 
fn=g n =0, and provide relations for evaluating the constants v n and A n . rj n thus can be ex- 
pressed in the form n n =v n {t—xlU n )+A n t , which when substituted into (5), gives as a general 
solution for $ mj $ n =xy{l + e n ) 2 -(v n IU n ) 2 -l]+h n (t) +K n . £ n is a solution of (9b) only when 
h(f) is a linear function of t. Thus g w can c onveniently be expressed as % n =— w n (t—xlU n ) + 
B n t+K n , where wJU n ='yl(l + € n ) 2 —(v n IU n ) 2 —l. The boundary conditions to be satisfied 
are £ p =0 at x=0, £ n =£„-i at x=U n t, and fc= -Wit+eiUit at x=Uit. These conditions 
require that the K n =0, and provide relations for evaluating the constants w n and B n . 

The solution for the transverse-wave region thus can be written; 
For the interval U 2 t<x<U 1 t, 

n^tkit-x/Ux), (27a) 



*i 



-Wnit-z/UJ-Wit + eJJtf, 



WitU^il + ^-fa/UJ*-!, 



€=€i 






da. 



For the interval U 3 t<x<U 2 t, 

v 2 =v 2 (t-x/U 2 )-t[v 1 (l-U 2 fU 1 )], 

i2 =-w 2 (t-x/UJ-tfa (l~U 2 fU 1 )+W 1 -e 1 U 1 ], 



W2/U 2 =J(l + e 2 )*-(v 2 /U 2 y-l, 



■e 2 . 



(27b) 

(27c) 
(27d) 

(27e) 
(27f) 

(28a) 
(28b) 
(28c) 
(28d) 
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For any interval U n+ it<x<U n t, 



Vn=v n {t-xlU n )+tit l v k (l-lJ w IU k ), (29a) 

k=l 

fc.= -w n (t-x/U P )-t V^w k {l-U k+l [U k )+W l -e l U^ (29b) 



v>n/U n =TKl + e n )'-(vn/U n )*- 1 , (29c) 



€=€ w . 

For the interval 0<x<U p t 1 



(29d) 



rj p =v p (t-x/U p )+t & Z) * v k (l-U k+1 /U k ), (30a) 



&=i 



£,= -w p (t-x/U p )-t ^^ 1 w k (l-U k+1 fU k )+W 1 -eJJ 1 ^ (30b) 

v> P /U p =j(l + e p )*-(v p fU p )*-l, (30c) 

€=€,. (30d) 

A point some distance down the filament thus is motionless until the elastic-wave front 
reaches it. Subsequently the strain at that point increases, and the point moves inward 
horizontally at increasing velocity until the tip of the transverse wave reaches it. The strain 
is then e t and the inward velocity W\. As the first segment of the transverse wave passes through 
the point the trajectory of the point turns upward in a direction inclined to the vertical. As 
subsequent elements of higher strain pass through the point the trajectory becomes less 
inclined, until at strain e p the direction of travel is vertical, parallel to the impact direction. 

The motion of filament material remains vertical in the region 0<x<U p t. There is no 
horizontal motion of material in this region because of the boundary-condition requirement 
at x=0 that £ p =0. In the region of the strained segments ei<€<e p , some flow of material 
takes place along the horizontal direction. It is this flow that causes the trajectory of a fila- 
ment point to be inclined from the vertical. The vertical and horizontal components of velocity 
of a filament point in this region are given by the coefficients of t in eq (29a) and (29b) . 

Although for generality each of the strained segments e n shown in figure 5 has been assumed 
to have a different angle 6 n of inclination to the horizontal, it will be proved that the angles B n 
are all equal and that the configuration of the transverse wave is a straight line. 

The displacements £ n , rj n , and the horizontal and vertical components of the tension are 
required to be continuous at each of the points x= U\t, . . .U n t, . . JJ P t. At the point x= U\t } 
fi and 77! are continuous by definition, and continuity of the horizontal and vertical tension 
components leads to eq (27d) and (27e). At x= U n t and at x= U p t, % n and t\ n are continuous by 
definition. Equating horizontal tension components gives 



if7,V (i+ " )i -te)- Mc '{«'- + 'ir»< i -% ! ) + «"- < ' t7 '] 



-ifeV ti+, --' )! -fe)'- M 4'"--' + 'SH i -% ! ) +,, ''-'' ff '} 

or with the help of eq (29c) 



[ifi;- Mf7 -]V (I+< - ) '-te) , -[ife- MC '']V (1 - | -'"- , "-feT- (31) 
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Equating vertical tension components gives 



l + e„Z7 : 



;-Mr,[,4;;(,-^)]= I ^^-w.[,, J+ '|'.,(.-|i)] 



or 



Manipulation of eq (31) and (32) gives the relation 

[-^-MJl] (l + 6j = [ T ^--MC7|](l +«,_,), (33) 



which reduces to 



_ / 1 (r.-r-,) , 



Equation (34) also holds when n=p. Dividing (32) by (33) gives a simple relationhip 

(l+e n )U n (1 + ^-JUn-t* 

This relationship also holds when n=p. Using (6a), eq (35) becomes 



(35) 



^7^-siniV (36) 



(l + e n )U n (l + e p )U p 
Applying (36) to (30c), (29b), and (29c), and using (6b), gives the relation 



At x=0 the boundary conditions rj p =Vt and % p =0 must be satisfied. The first of these 
conditions gives 

V=v p +~f± £ (U k -U k ^)=v p + n . V \ TJ k= J2 \l+e k )(U k -U k+1 ) 

k=i Uk (l-tepjUp k=l 

= nJ\ TT I *£ 1 ^ +1 (€* +1 -€,) + (l+€ 1 )?7il (38) 
(l-f-e p )U p L fe=i J 



By the use of (38) , eq (36) can be expressed as 

». v 



(l + * n )U n 2+(l+ ei )Z7! 
where 2 is an abbreviation for 



= -sin0 B , (39) 



2= l]V« +1 (ft +I -ft), (40) 

Similarly, eq (37) becomes 

i^O+^vHSwiWH''-- (4i) 

Equation (39) is equivalent to a statement that each sin 6 n has the same value given by 

sin e »=- s+( i+ fl )c/r sin * ( 4 ' 2 ) 
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Similarly, eq (41) states that each cos 6 n has the same value given by 

co8g„= A /l- [ s+(1 ^ <i)£/ J =cos 6. (43) 

The configuration of the transverse wave, therefore, is a straight line. 
The boundary condition J„=0 requires that 

w P + g ^(l~^± i )+Fi--6 1 t7 1 =0. (44) 

It can be shown that 

, ^ t IH 1 -^)=ii;( ] +M fc S Vi+l(ft+1 - fe)+(1+£l)cri ]- C71 

= [S+(l+€ 1 )C7 I ]cos^-C/ 1 



w„ 



=V[2 + (l+€ 1 )i!7 1 ] 2 -F 2 -Z7 1) (45) 

so that eq (44) can be put into the form 



y=V[2+(l+a 1 )t/ 1 ] 2 -[(l+ ei )t7 1 -W 1 ] 2 . (46) 

Equations (46) and (42) \;an be combined to obtain two different expressions for sin 9 and cos 



smg =-V nW+e^J ' (47) 

cose ' s+a+nWi ^' 

Equations (31) through (48) can now be used to modify the solutions given earlier. Thus 
(30a) and (30b) become 

-K ( -w¥fe}- (49a) 

= { (i+6,)cos e-i }x=- { 1 -( 1 +^)[ s+(i+ 1 ^ 1 ] } a; - (49b) 

and (29a) and (29b) become 

„=<{ 8 ,+*gv«(i-^)}-^ 

= '{(H^[¥ P -' < -«- a) +' 1 +' JP -]}- Z+(t+ ) .,V 

1 1'"' s+ t i+«x7. -J'-Ls+o+^J'J- 



(50a) 
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^=-^W n + k= ± 1 Wk (^l-^^+W i -e 1 uXt+^X 
=- /PgVt+.te+i-ftRa + eOtfiJcos 6-l(l + ei )U 1 -W l ]\t-{l-(l + e n ) cos 6}x 

The equations just given are for the case in which the tension-strain curve is approximated 
in the region e x <e<e p by a finite number of straight-line segments. The solution sought is 
the limit solution when the number of straight-line segments approaches infinity and the 
length of each segment approaches zero. This limit solution is given in the following section. 

7. Tabulation and Interpretation of the General Solution 
For the interval 0<x<U p t, 

"= r {>-sW+Su}- (510) 

e=e p . (51c) 

For the interval U p t<x<U 1 t } 

For the region Uit<x<^ oo the solution is given by eq (25) and (26). The symbols used in the 
above equations have been defined previously, except for 



s -i' t 'yfM da = w ~ w ^ (53) 

and 

«W- [\lh d ld-=^)-W, (54) 



The general solution given in this section differs from the simple solution given previously 
in that the strain in the transverse wave need not be constant. At low impact speeds the 
equations of the general solution reduce to those of the simple solution. The transverse wave 
then propagates in a region of constant strain in the wake of the plastic-wave front. How- 
ever, for a sufficiently high impact velocity and a tension-strain curve of the right shape, the 
case in which the transverse-wave front and the plastic-wave front coincide occurs. The 
velocity of propagation of these two wave fronts is TJ\ and the maximum strain in the filament 
is cj. If the tension-strain curve is known, the values of JJ\ and e x for this special case can be 
found from eq (27e) . 
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At higher transverse impact-velocities the plastic-wave front propagates more slowly than 
the transverse-wave front. A region of variable straining therefore occurs in the transverse 
wave, and the general solution must be used. The shape of the transverse wave remains a 
straight line, but the propagation velocity no longer increases with increasing impact velocity. 
It remains constant at the value Z7i. The velocity of propagation of the maximum strained 
portion of the transverse wave U p is the same as the plastic-wave velocity C p and is computed 
from eq (2). In order to find the value of the maximum strain, or constant strain e p in the wake 
of the plastic-wave front, the relation 



V=4S+(l+€ 1 )U 1 ?-[(l+e l )U l -W l ] 2 (55) 

can be used. The value of S depends upon e p , so that if the tension-strain curve is known, the 
value of € p needed to give a required V can be found by a process of trial and error. 

At the transverse-wave front, material of the filament is flowing inward horizontally with 
velocity W\. At the point in the transverse wave where the strain attains its maximum value 
e p , material is flowing inward along the filament direction with velocity W. The quantity S 
represents the increase in velocity over W x due to the strain increase in the transverse-wave 
region. 

The motion of filament points in the wake of the plastic-wave front is most easily visualized 
as follows; Suppose there is a second filament, inextensible and fixed to the impacted filament 
at a point a great distance in advance of the elastic-wave front. Let the end of this filament 
coincide with the impact point. After impact let the filament be constrained to move alongside 
the impacted filament but without extending so that its end no longer coincides with the im- 
pact point, but recedes from it. To an observer at this filament end the impacted filament 
would appear to have motion with velocity W. 

Although flow of filament material in the wake of the plastic wave occurs, there is no 
horizontal component of the flow velocity. A horizontal component of the flow does occur 
in the region where the strain varies: i. e., where ei^ef^. The magnitude of this velocity 
component is much smaller than Wi, the horizontal flow velocity in advance of the transverse 
wave, and decreases to zero where e equals e p . 

At low impact speeds, where the simple theory applies, the angle that the transverse 
wave makes with the horizontal can be found from the equation 

As the impact velocity is increased, cos decreases to the value 

.a+eoc/i-w; 



COS $i = 



(l+«l)tfl 



(0 is a negative angle and attains a greater negative value as the impact velocity is increased.) 
At higher impact speeds, where only the general theory applies, 6 can be found from the 
relation 

cozC tt+'i)^-^ (57) 

The values of €i, Z7i, and W\ are constants, but S increases as the impact velocity is increased. 
Cos continues to decrease, and the value of becomes greater negatively. Thus, although 
the progagation velocity of the transverse wave is constant for all impact speeds above a certain 
value, the slope of the transverse wave continues to change. 

Figure 6, a, shows the configuration of the filament after a transverse impact of sufficiently 
slow speed that the simple theory applies. Figure 6, b, shows the configuration after a higher 
speed impact where only the general theory applies. As pointed out previously, the formula 
for U expresses the transverse wave velocity relative to a Lagrangian coordinate system. 
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Figure 6. Configuration of filament t seconds offer transverse 
impact. 

a. Simple theory; b. Genera] theory, 



i— [{!+€,) U,-WQt - 



In laboratory coordinates the distance that the transverse wave travels along the filament in 
unit time is (l + e p )?7 for the simple theory, and S-\- (l + €i)i7i for the general theory. The cor- 
responding distances that the wave front travels in the horizontal direction are (l J r€ p )U—W 
and (l + €i)C7i — W 1} respectively. These distances are labeled in figure 6, a and b. Thus it is 
seen how some of the relationships derived above could have been obtained very simply from 
a knowledge of geometry and an intuitive understanding of the problem. 

8. Examples of the Theory 

In order to illustrate the theory, the behavior of a filament after transverse impact is 
computed for two different hypothetical materials. The first material obeys Hooke's law. 
Its tension-strain curve is given by 

T 



M~ 



--C 2 



(58) 



As this tension-strain curve is a straight line of constant slope, the simple theory applies. 
The elastic- and plastic-wave fronts both propagate with velocity C. U, V, and W are given 

by 



u=c^ 



l + e p 



V=C^ p {l + e v )-Ut v (l + e p )-e p ]\ 



W=Ce p . 



(59) 

(60) 
(61) 



Curves of U/C, V/C, and W/C as functions of e p for the Hookean material are plotted in figure 7. 
The solution of a transverse-impact problem may be found from these curves by reading off the 
values or e p , U/C, and W/C corresponding to a particular V/C. The value of C is determined 
from the Young's modulus and density of the unstrained material. 

In order to illustrate a case in which the general theory applies, computations were made 
for a material whose tension-strain curve is given by 3 



j£= 14 tanh 20e- 



(62) 



For this curve, plotted in figure 8, the units of T/M are arbitrary, and the strain e is expressed as 
a fraction of unity. The slope of the curve, initially 280 arbitrary units per unit strain, con- 
tinually decreases and approaches a value of zero for very large strains. The slope d{TjM)jd€ 
and the quantity T/M(l+e) become equal to 12.308 at a strain of 0.11221. When e is 0.20, 
the highest strain for which computations were made, the slope is 0.37547. In the region 
0<a;<0. 11221 the simple theory applies. In the region 0.11221 <e<oo the general theory 
applies. 

The elastic-wave velocity is C e =^j2SQ= 16.733 arbitrary velocity units. The plastic- 
wave velocity is given by 



0p = v i(S) = ^ 280 sech2 2 ° €p= ° e sech 20e 



(63) 



3 The functional relationship (62) was chosen because tables were available that simplified the calculations. The numbers 14 and 20 have no 
significance, and were chosen only for convenience in plotting the curve. Other relationships involving more familiar functions could have been 
used^such as, T/M=a+b€— ■y/a 2 +b 2 e 2 , where a and b are arbitrary constants. 
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Figure 8. Plot of tension-strain curve given by 
(T/M) = 14 tanh20e. 

Point of tangency with dotted line indicates strain and tension at which the 
slope of the curve equals T/M(l+c). 



Figure 7. Curves of U/C, V/C, and W/C as functions of e p for 
a filament having a linear tension-strain curve. 

In the region 0<e p <0. 11221 the transverse-wave velocity is given by 



TT _ I T _ n / 0.05 tanh 20e p 

= VM(l + € ,rS/ 1 + e, 



(64) 



When e p is equal to or greater than 0.11221, the transverse-wave velocity is constant at the 
value £7 1 =0.20967ft. 

Filament material in the wake of the plastic-wave front flows inward at a velocity given by 



W= | €V Jji^ da = omC e tan-^sinh 20«„)- 



(65) 



Figure 9 gives the curves of U/C e , V/C e , and W/C e as functions of e p for this material. 
The dotted lines give portions of the corresponding curves for the Hookean material. 

Other behavior features of a material having the nonlinear tension-strain curve of eq (62) 
are given in figures 10, 11, 12, and 13. Figure 10 shows the strain distribution along the fila- 
ment after transverse impact at various speeds up to 0.177ft. The dimensionless parameter 
xJC e t is used as the abscissa. A particular value of this parameter, 0.6 for instance, can be 
interpreted as a distance down the fiber 0.6 as far as the elastic wave has propagated. The 
dotted line gives the strain at the transverse-wave front as a function of the location of the 
transverse-wave front. Thus the strain in the filament after a transverse impact at velocity 
of 0.151ft increases as shown up to a strain of 8 percent between the elastic- and plastic-wave 
fronts, and then remains constant. The plastic-wave front propagates at a velocity of C p = 
0.388(7*, and the transverse-wave front at a velocity of U= 0.207ft for this case. 

Figure 11 shows the horizontal component of the inward-flow velocity for different points 
along the filament. The location of a point on the filament is designated by the same distance 
parameter used in figure 10. Consider, for instance, the case in which the filament is impacted 
at velocity 0.173ft. The strain in the wake of the plastic-wave front is 16 percent. The 
flow velocity at various points along the filament between the elastic-wave front and the 
transverse-wave front at #/ft£=0.210 increases to the value 0.0680ft as shown. A& the fila- 
ment is horizontal in this region all of the flow is horizontal and directed inward. In the wake 
of the transverse-wave front the filament has changed its direction. The horizontal component 
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Figure 10. Strain distribution along filament having tension- 
strain curve given by (T/M) = 14 tank 20 e, after transverse 
impact. 

If the filament is impacted, for example, with velocity 0.151C, the strain In- 
creases to a value of 8 percent in the region between the elastic- and plastic-wave 
fronts, and then remains constant. The plastic-wave front in this case propagates 
at a velocity of C P =0.388C«, its position along the filament in distance parameter 
units is thus 0.388. The position of the dotted line gives the position of the trans- 
verse-wave front, which in this case propagates at velocity (7=0.207 C„ 
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Figure 9. Curves of U/C e , V/C e , and W/C e as functions 
of e v for a filament with tension- strain curve given by 
(T/M)=14tanh20e. 

Dotted lines show portions of similar curves- for filament having linear stress- ° "' ' 2 ,3 «.-«j;l- b ,bLtm 8 « /^ J 

strain curve. distance parameter, x/Ce' 

Figure 11. Horizontal component of flow velocity at various 
points along the length of filament having tension- strain curve 
given by (T/M) =14 tanh 20e, after transverse impact. 

A filament impacted at velocity V=0.173C e , for example, has a strain in the 
wake of the plastic-wave front of 16 percent. The inward-flow velocity at various 
points along the filament between the elastic-wave front and the transverse-wave 
front at z/Cy= 0.210 increases to the value 0.0680C e . In the wake of the transverse- 
wave front the filament has changed direction, and the horizontal components 
flow velocity, therefore, decreases abruptly to the value 0.00456%, and then de- 
creases gradually to zero at x/C e t=Q.Q8l, the location of the plastic-wave front. 

of the flow velocity, therefore, decreases abruptly to the value 0.0045CV It then decreases 
gradually to zero at xj C^= 0.081, the location of the plastic-wave front. Between the plastic- 
wave front and the impact point the horizontal component of the flow velocity is zero. 

Figure 12 shows how the angle of the transverse wave varies with the strain e p in the 
wake of the plastic-wave front. The solid line gives calculated values of cos and for the 
material with the nonlinear tension-strain curve. The dotted line gives the values calculated 
for the Hookean material. The angle increases negatively very rapidly for small values of 
€ p . Thus even at low values of e p and correspondingly low impact velocities, a well-defined 
tent-shaped wave may be expected. (For instance, if C e = 2,400 m/sec and F=0.026<7 e =62.5 
ra/sec, e p =0.5 percent and 0=-— 21.75°.) At higher values of e v the configuration of the 
transverse wave does not change much. This is illustrated in figure 13, where several con- 
figurations after impact are plotted for the filament having the nonlinear stress-strain curve. 
The locations of the plastic-wave front on several of these configurations where e p exceeds 
0.11221 are given by the black dots. 
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Figure 12. Angle which transverse wave makes with hori- 
zontal plotted as a function of the strain (e p or ey) at the trans- 
verse-wave front. 

r Solidjline: Filament having tension-strain curve given by (T/M)=U tanh 20e. 
Dotted line: Filament having linear tension-strain curve. 




Figure 13. Configurations after transverse 
impact for filament having tension-strain 
curve given by (T/M)=14 tanh 20e. 

a. Before impact. 

b. e p =0.04, F=0.110 C e , 0=-36.4°. 

c. e p = .08, V= .151 C«, 0=-42.5°. 

d. e p = .12, V= .166 C e , 0=-45.3°. 

e. € P = .16, V= .173 C e , 0=-46.4°. 

f. e p = .20, V= .177 C e , 0=-46.9°. 

Black dots in last three configurations show positions 
of plastic-wave fronts. 
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